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ABSTRACT 

We present multi-line and continuum observations of the circumstellar envi- 
ronment within 10 4 AU of a sample of protostars to investigate how the effects of 
outflows on their immediate environment changes over time. 12 CO (1-0) emission 
probes the high-velocity molecular outflows near the protostars and demonstrate 
that the outflow opening angle widens as the nascent star evolves. Maps of the 
13 CO (1-0) and HCO + (1-0) outflow emission show that protostellar winds erode 
the circumstellar envelope through the entrainment of the outer envelope gas. 
The spatial and velocity distribution of the dense circumstellar envelope, as well 
as its mass, is traced by the C 18 O(l-0) emission and also displays evolutionary 
changes. We show that outflows are largely responsible for these changes, and 
propose an empirical model for the evolution of outflow-envelope interactions. 
In addition, some of the outflows in our sample appear to affect the chemical 
composition of the surrounding environment, enhancing the HCO + abundance. 
Overall, our results confirm that outflows play a major role in the star formation 
process through their strong physical and chemical impacts on the environments 
of the young protostars. 

Subject headings: stars: formation — ISM: jets and outflows — ISM: evolution 
— stars: pre-main-sequence — ISM: molecules 
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1. Introduction 

Young stellar objects ranging in age from approximately 10 3 yr to more than 10 6 yr are 
known to power outflows that entrain the surrounding gas of the parent molecular cloud. 
Protostellar winds originate close to the surface of the forming star and must interact with 
the dense envelope surrounding the protostar-disk system (e.g., Konigl & Pudritz 2000; Shu 
et al. 2000). These circumstellar envelopes, with sizes of 10 3 to 10 4 AU, are the primary mass 
reservoirs for young stellar objects (YSOs). As a result, perturbations caused when outflow 
momentum and energy are injected into the protostellar surroundings can have a major 
impact on the mass- assembly process. In fact, outflows are most likely to be responsible 
for the clearing of material from the protostellar surroundings, a process that could result 
in the termination of the infall phase (e.g., Velusamy & Langer 1998) and affect the star- 
formation efficiency in the cloud (e.g., Matzer & McKee 2000). However, there is still no 
general understanding of how the outflow's effect on its immediate surroundings changes as 
the protostar evolves, as studies of the outflow-envelope interaction have been limited to a 
few specific cases. 

In the nearest star-forming regions, circumstellar envelopes around low-mass young 
stellar objects have angular sizes of about 10" to 30". A detailed study of the outflow 
impact on the envelope requires observations that distinguish between the various velocity 
patterns that can exist in different density regimes in these envelopes (i.e., infall, outflow, 
rotation, and turbulence). Observations at high angular and velocity resolution of a variety of 
molecular species that trace a range of density and kinematic regimes are therefore essential. 
Also, in order to examime how the outflow-envelope interaction evolves, observations of the 
environment around protostars of different ages are needed. 

Here, we present the results of a systematic survey of the circumstellar environments 
of a sample of YSOs with outflows. Multi-line and continuum high spatial and spectral 
resolution observations from the Owens Valley Radio Observatory (OVRO) millimeter-wave 
array enable investigation of the outflows, circumstellar gas and evolution of the outflow- 
envelope interactions, on scales of ~ 10 3 AU. Considerably more detailed studies of two of 
the sources in our sample, L 1228 and RN0 43, have already been published (Arce & Sargent 
2004; 2005, hereafter Paper I and Paper II). 

2. The Sample 

Our observed sample consists of nine low-mass protostars at different evolutionary 
stages, and include: 1) three class sources — deeply embedded protostars in the earli- 
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est formation stage, with ages between ~ 10 3 and ~ 10 4 yr (see Andre et al. 2000); 2) three 
class I sources — protostellar objects still embedded in their parent cores but older than the 
class sources, with ages ~ 10 5 yr (see Lada 1999); and 3) three class II sources, with ages 
of a few ~ 10 6 yr (Lada 1999). All are within 500 pc of the Sun and are known outflow 
sources based on either millimeter or optical observations. Approximate protostellar ages 
are obtained from the bolometric temperature, Tboi, the temperature of a blackbody having 
the same mean frequency as the observed spectral energy distribution (SED) of the YSO 
(see Myers & Ladd 1993; Chen et al. 1995; Ladd et al. 1998). Sources in our sample have 
T bo i ranging from 40 to 2700 K, corresponding to ages from ~ 10 3 to ~ 10 7 yr. Table 1 lists 
the source names, coordinates of the observed fields, and other relevant parameters. A brief 
description for each source is given in the Appendix. 



3. Observations 

The observations were carried out using the OVRO millimeter array of six 10.4 m 
telescopes, during three observing seasons between 2001 September and 2004 May. For all 
sources in our sample, a single pointing (i.e., no mosaicing) and two different correlation 
configurations were used to provide images on a total of six different molecular lines and two 
4 GHz-wide continuum emission bands. The area on the sky covered by each map is the 
OVRO primary beam, 65" x(100/V) (FWHM), where v is the observed frequency in GHz. 

The first digital correlator configuration provided simultaneous measurements of 12 CO(l- 
0), 13 CO(1-0), and C 18 O(l-0) at 115.27, 110.20, and 109.78 GHz, respectively, and 4 GHz 
bandwidth continuum emission at 2.7 mm. With the local oscillator (LO) frequency set 
to 112.736 GHz, the 12 CO line was observed in the upper sideband with a channel width 
of 0.325 km s _1 over a 20 km s" 1 bandwidth. In the lower sideband, the 13 CO and C 18 
lines were each observed at a velocity resolution of 0.17 km s" 1 over a 5 km s _1 bandwidth. 
Observations were obtained in two to five different array configurations depending on the 
source and with baselines extending from 15 m to a maximum of 480 m. 

Simultaneous observations of HNC(l-O) F=2-l, HCO+(1-0), and H 13 CO+(1-0) at 90.66, 
89.19, and 86.75 GHz, respectively, and of 4 GHz bandwidth continuum emission centered 
at 3.4 mm were obtained with the LO frequency tuned to 87.971 GHz, so that the HNC 
and HCO + lines were in the upper sideband and H 13 CO + in the lower. All three lines 
were observed with two or three different array configurations with spectral resolutions of 
~ 0.1 km s^ 1 over bandwidths of 6 km s _1 , and baselines from 15 m to a maximum of 
120 m. Only maps of the HCO + emission are discussed here, as they are relevant to the 
present discussion. 
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A bright (> 1 Jy) quasar near each observed source was used for phase and gain calibra- 
tion. Flux calibration relied on observations of 3C273 and Uranus, with resulting uncertain- 
ties of about 20 to 30%. Visibility data was edited and calibrated using the MMA software 
package (Scoville et al. 1993) and images were produced using the MIRIAD software package 
(Sault et al. 1995) and its CLEAN algorithm, with "robust" weighting of the visibilities. For 
each map the robust parameter was chosen to optimize both spatial resolution and sensitiv- 
ity to extended structures, as well as the signal-to-noise ratio. Table 2 provides a list of the 
relevant parameters for the molecular line and continuum emission data. 

4. Results 

The gas surrounding forming stars is generally characterized by complex kinematics and 
steep density gradients and observations of different molecular lines are necessary to trace 
different density and kinematic regimes. We use the optically thin 3 mm dust continuum 
emission to pinpoint the position of the sources in our sample and to study the inner circum- 
stellar environment. The 12 CO(1-0) emission probes the low-density, high- velocity outflow 
gas, while the less optically thick 13 CO (1-0) reflects the behavior of the low- velocity, slightly 
denser outflow, as well as the impact of the protostellar outflow on the outer regions (> 2000 
AU) of the circumstellar envelope. Observations of the C 18 O(l-0) line trace the inner dense 
parts of the envelope surrounding the protostar and in some cases may be used to probe 
the effect of the outflow on the high-density gas. We originally observed HCO + (1-0) with 
the expectation that this would trace infall motions of the gas surrounding the young stars. 
Since HCO + (1-0) kinematics appear to be dominated by outflow motions, we have used these 
results to examine the physical and chemical impact of the outflow on the surrounding gas. 

4.1. Continuum Emission 

Millimeter continuum thermal dust emission at 2.7 and 3.4 mm was detected towards 
all the sources in our sample. Maps of this continuum emission were constructed so that 
both bands had similar synthesized beams (see Table 2). Since the linear resolution of our 
continuum observations ranges from about 700 to 2700 AU, it is very probable that both the 
emission from the envelope around the protostar and that from an unresolved circumstellar 
disk contribute to the total detected emission. We expect that in the class sources more 
than 80% of the emission comes from the inner (< 2000 AU) circumstellar envelope. For 
the class I and II sources, we expect the contribution from the disk to the total detected 
emission to rise from about 30% to 90% as the source evolves (e.g., Hogerheijde et al. 1997a; 



- 5 - 



Looney et al. 2000). 

We estimate the total mass of the gas and dust associated with the continuum emission, 
M cont , from F V D 2 / B u (T d )K u (e.g., Beckwith et al. 1986). Here, F v is the continuum flux 
density at frequency v, D the distance to the source, B v the blackbody (Planck) function, T d 
the dust temperature, and k v is the dust mass opacity coefficient, with the usual form, k v = 
0.1(z//10 12 Hz)' 3 cm 2 g" 1 (e.g., Beckwith & Sargent 1991). The value of (3 was determined by 
fitting the measured flux densities at wavelengths between 0.62 and 3.4 mm with a power law, 
F v oc is( 2+/3 \ using measurements from our observations and the literature (see references in 
Table 1). For sources with less than three available flux measurements, we adopt (3 = 1. The 
dust temperature was estimated from greybody fits to the spectral energy distribution of the 
source (e.g., Chini et al. 1997) at wavelengths greater than 60yum and using the estimated 
(or assumed) value of (3. The measured fluxes at 2.7 and 3.4 mm, F 2 . 7 and F 3A , respectively, 
and the estimated values of (3, T d) and M cont are listed in Table 3. 



4.2. 12 CO(1-0) outflows 

In Figure 1 we present integrated intensity maps of the 12 CO(1-0) emission associated 
with the observed sources. The central envelope velocity, v env , and the velocity ranges of 
integration are given in Tables 8 and 4, respectively. Blueshifted and redshifted components 
are detected for all sources; velocities deviate by at least 0.5 km s _1 , and usually by more 
than 1 km s" 1 , from v env . Since our interferometer observations filter out structures on scales 
greater than ~ 30", we are confident that the emission seen in Figure 1 traces molecular gas 
accelerated by the opposing outflow lobes close to the protostellar sources, and not ambi- 
ent material. However, filtering of large scale structures in our interferometer observations 
hampers our ability to fully recover the flux from the outflow gas present in the observed 
area. 

The amount of filtered or missing flux in our OVRO data depends on the observed 
molecular outflow's structure; a smaller fraction of the total flux will be recovered from 
molecular outflows composed of smooth large-scale structures, compared to the fraction of 
the total flux recovered from molecular outflows that are mainly made of compact structures. 
For three sources in our sample, we compared published 12 CO data from the Five College 
Radio Astronomy (FCRAO) 14 m telescope (Tafalla & Myers 1997; Lee et al. 2000; Arce & 
Goodman 2001) with our OVRO data averaged over a 46" x 46" area, the approximate beam 
size of the FCRAO 12 CO observations. We find that for L 1228 our OVRO observations 
recover more than 60% of the flux at all outflow velocities, and more than 90% of the flux at 
high outflow velocities. On the other hand, in HH 300 and RNO 91 we recover 15 and 20%, 
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respectively, of the single-dish flux at outflow velocities. Evidently, the amount of missing 
flux is different for each molecular outflow. Here we do not attempt to correct for the missing 
flux in our OVRO data, as we lack the necessary single-dish data to do so for all sources. 
We stress that regardless of how much flux is missing in each of the molecular outflow maps 
in our sample, our data are sensitive to the small-scale outflow structures — important for 
the study of the outflow-envelope interactions — and our maps reliably show the outflow 
morphology close to the source. 

4.3. 13 CO(1-0) outflows 

In Figure 2, integrated intensity maps of the 13 CO emission towards our sample sources 
show gas that has been accelerated by the protostellar outflow and is slower and more dense 
than the 12 CO features. Like the 12 CO features, the morphology of the 13 CO blueshifted and 
redshifted outflow lobes varies substantially among sources. Nevertheless, with the exception 
of the red 13 CO lobe of GK/GI Tau, the locations and positional offsets of the 13 CO and 12 CO 
outflow lobes correspond well. In several sources the 13 CO outflow emission is concentrated at 
the outer edges of the 12 CO outflows, and appears to delineate the dense outflow cavity walls. 
The blue lobes of L 1228 and RNO 129 are particularly nice examples of this behavior but 
it is also suggested for RNO 43. In other sources the 13 CO outflow concentrates at positions 
slightly offset from the protostar in the same direction as the corresponding 12 CO outflow 
lobes, and generally have opening angles larger than that seen in 12 CO. This morphology is 
best illustrated by HH114mms, HH300, the red lobe of RNO 129, and the blue lobes of T 
Tau and IRAS 3282. Close to these sources, blueshifted and redshifted 13 CO traces dense 
gas that has been propelled by the protostellar outflow, as well as entrained dense gas in 
the outer outflow walls. On ocassion, 13 CO is concentrated outside the limits of the 12 CO 
outflow, as in the red lobe of T Tau and the blue lobe of GK/GI Tau, perhaps due to the 
presence of a dense shell produced by the interaction between the protostellar outflow and 
the surrounding environment, and similar to the larger scale morphology observed in other 
sources (e.g., Arce & Goodman 2002a). 

Similar to the 12 CO data, we expect that our 13 CO interferometer maps filter a sub- 
stantial amount of flux (see §4.2). For example, comparing single-dish (FCRAO) data of 
L 1228 and HH300 (Tafalla & Myers 1997; Arce & Goodman 2001) with our OVRO maps, 
we find that the interferometer 13 CO observations recover only 10% and 15%, respectively, 
of the flux in the outflow wings. This suggests that our OVRO observations filter most of 
the 13 CO (large-scale) cloud emission and the interferometer maps mainly show the (small- 
scale) emission associated with the envelope and the regions where it is being impacted by 
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the underlying protostellar wind. 

4.4. Molecular Outflow Parameters 

The 12 CO line is expected to be optically thick at velocities close to the ambient envelope 
velocity and optically thin at high outflow velocities. Thus, we use 12 CO emission to estimate 
the mass of the low-density, high-velocity component of the molecular outflow and 13 CO 
emission to account for higher density, lower velocity material (see papers Papers I and II 
for more details). At outflow velocities where the 12 CO line is extremely optically thick and 
there is detectable resolved 13 CO outflow emission, we derive the outflow mass from the 
13 CO emission, as it is is expected to be optically thin at these relatively high velocities far 
from the line core. Outside the velocity regimes where we use 13 CO, we estimate the outflow 
mass from the 12 CO emission. 

We estimate the outflow mass, M out , from the 12 CO and 13 CO data assuming LTE 
and a gas temperature of 25 K. We adopt [H 2 ]/[ 13 CO] = 7 x 10 5 (Frerking et al. 1982) 
and [ 13 CO]/[ 12 CO] = 62 (Langer & Penzias 1993). We obtain the outflow momentum, 
P out = T, v M out (v)v, and kinetic energy, E out = ^Y, v M out (v)v 2 , for each source and list their 
values in Table 5. In addition, we calculate the mass outflow rate, M out = T, v M out (v )/r, and 
the flow momentum rate, Fqo = ^vM out (v)v/r (see Table 6), where r = v/r is the dynamic 
time scale of the outflow gas, v is the outflow velocity, and r is the distance from the source 
to the outflow gas. In Table 6 we also show the dense gas mass outflow rate, M^ ense = 
T, v Mdense(v)/r, where M& ense is the outflow mass traced solely by the 13 CO emission. 

For all mass-related outflow parameters listed in Tables 5 and 6 (except for M dense , where 
we only use the 13 CO emission) the final estimate is the summation of the estimates from 
both the 12 CO and 13 CO emission. These values should be treated as lower limits, since no 
correction for missing flux, opacity or for outflow inclination angle have been applied. Care 
should also be taken when comparing M out , P ou t, and E out for the individual outflows since 
the areas observed are in effect, the linear sizes of the interferometer primary beam, and 
differ according to source distance. Nevertheless, to the first order M out , F C o, and M dense 
do not depend on the observed area (Henning & Launhardt 1998). 

4.5. HCO+ outflows 

In Figure 3 the HCO + emission at | v env |> 0.5 km s -1 appears to coincide with the 
molecular outflow lobes. Indeed, outflowing HCO + (1-0) is observed in all sources except 
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GK/GI Tau, although not all show both blueshifted and redshifted emission. In general, 
this HCO + is confined to the outflow cavity walls. In this respect it echoes the behavior of 
13 CO but, as can be seen in Table 4, the velocity range exceeds that of the 13 CO outflow 
emission in most sources, and, as can be seen in Figure 3, in some sources the HCO + outflow 
emission does not trace the full extent of the outflow cavity walls. As in the 13 CO data, 
we expect that our interferometer observations have filtered much of the extended emission, 
and that the HCO + outflow maps show the compact structures in the envelope where the 
gas has been influenced by the prostellar outflow. 

The HCO + abundance in the outflow, Xhco+ > should indicate whether outflow-triggered 
chemical processes have affected the surrounding gas. The 13 CO outflow emission along the 
outflow cavity walls and offset from the source position is probably optically thin, since 
there is no C 18 at these velocities and positions. Thus, it provides a reliable estimate of 
the molecular hydrogen column density, Nh 2 , using the same assumptions as in § 4.4. The 
HCO + column density, N HC o+, can be obtained assuming the emission is optically thin, and 
leads to the outflow HCO + abundance using Xhco+ — ^hco+ I 'Nh 2 , where Nh 2 (from the 
13 CO data) and N HC o+ are obtained over the exact same area and velocity range. The uv- 
coverage and beam sizes of our 13 CO and HCO + data are very similar. Thus, these maps are 
sensitive to similar structures in the emission and may be used to obtain a reliable estimate 
of the HCO + abundance in the outflow and outflow-envelope interface, close to the protostar. 
Given that an estimate of X HC o+ c & n om Y be obtained in regions where the HCO + and 13 CO 
emission coincide in position and velocity, and that in some sources the HCO + emission aries 
from positions and velocities where 13 CO is not detected, we can obtain only a lower limit 
to X HC o+ for most of our sources (see Table 7). 



4.6. C 18 O(l-0) envelopes 

Integrated intensity maps of the C 18 emission associated with our sample are shown 
in Figure 4. For the class and class I sources, C 18 emission usually coincides with 
continuum peak positions. Elliptical gaussian fits indicate geometric mean sizes (after beam 
deconvolution) of about 2000 AU (FWHM) independent of evolutionary stage. An exception 
is RN043 for which the C 18 emission region size is about 5600 AU (FWHM) (see Table 8). 
The compact morphology and spatial coincidence with the dust continuum peaks strongly 
suggest the C 18 emission traces the dense circumstellar envelope gas in these class and 

I sources. By contrast, C 18 emission is significantly offset from the positions of our class 

II sources, and may trace the remnants of the high-density (n > 10 4 cm -3 ) circumstellar 
envelope gas that has been mostly cleared by the protostellar outflow (see §5.4). 
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Comparing the total C 18 flux obtained from single-dish (FCRAO) data of two sources 
(HH300: Arce & Goodman 2001; IRAS 3282: Hatchell et al. 2005), with the flux recovered 
by our interferometer data (smoothed to the same beam size), we estimate that the OVRO 
observations trace 1 to 4% of the total flux traced by the FCRAO observations. In addition, 
the velocity distribution of the emission in the single-dish data is very different from that 
of our OVRO data. These properties support the idea that the two sets of data trace 
entirely different structures: the C 18 emission detected in our OVRO data arises from the 
compact structure of dense gas surrounding the protostar (i.e., the circumstellar envelope, 
or its remnants), while the single-dish C 18 observations trace the surrounding (large-scale) 
cloud. Thus, our C 18 OVRO data can be used to obtain an estimate of the circumstellar 
envelope mass (e.g., Hogerheijde et al. 1998). 

We estimate the masses of the dense circumstellar envelope, M env , assuming optical 
thin C 18 emission, LTE, and an abundance ratio of [C 18 0]/[H 2 ] = 1.7 x 10~ 7 (Frerking et 
al. 1982). We adopt an excitation temperature, T ex , of 25 K, since T ex probably lies between 
10 to 40 K; the lower limit is given by the "typical" temperature of an isolated dense core 
(Jijina et al. 1999), while the upper is the "typical" Td of our sources. Envelope masses and 
other properties are listed in Table 8. 

5. Discussion 

5.1. Outflow Opening Angle 

The CO outflows displayed in Figure 1 have very different morphologies, due in part 
to their various surrounding environments. Nevertheless, there are significant trends. The 
molecular outflows from the youngest protostellar sources (class 0) tend to be highly colli- 
mated with either jet-like morphologies or cone-shaped lobes with relatively small opening 
angles of less than 55°. As discussed in Paper II, the morphology and kinematics of one of 
these, the RNO 43 outflow, strongly suggest that it results from the entrainment of ambient 
gas by an underlying collimated jet. Likewise, the significant collimation of the HH 114mms 
molecular outflow, as well as the fact that the outflow velocity increases with distance from 
the source, indicate that it is also driven by a jet or by a very collimated, momentum- 
conserving, wind (e.g., Lee et al. 2001). The cone-shaped morphology of the IRAS 3282 out- 
flow resembles that observed in other class outflows such as L 1448 (Bachiller et al. 1995), 
L1157 (Gueth et al. 1996), and HH211 (Gueth & Guilloteau 1999). These cone-shaped 
morphologies have been successfully fitted with jet bow shock entrainment models (e.g., 
Raga & Cabrit 1993). By analogy, we propose that the outflow lobes in IRAS 3282 represent 
the walls of a cavity produced by a jet bow shock. The limited velocity range of our 12 CO 
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observations (see Table 4) precludes our detecting the high velocity collimated component 
of IRAS 3282 that has been observed in single dish maps at outflow velocities of ~ 50 km 
s- 1 by Bachiller et al. (1991). 

The sightly more evolved (class I) sources in our sample have 12 CO outflows with lobe 
opening angles of more than 75°. This wide angle morphology is also seen in 13 CO (Figure 2) 
with opening angles for L 1228, HH 300, and RNO 129, which are as wide as, if not wider than, 
their 12 CO counterparts. Near-infrared and optical scattered light nebulae associated with 
these sources (Bally et al. 1995; Reipurth et al. 2000; Movsessian & Magakian 2004) display 
morphologies that are very similar to those of the molecular outflow lobes, consistent with 
the high-velocity 13 CO emission tracing the walls of cavities produced by the protostellar 
outflow. 

By contrast, outflows from the most evolved young stars in our sample, the class II 
sources, have no definite shape or structure. Blueshifted and redshifted gas within about 
6000 AU of RNO 91, T Tau, and GK/GI Tau shows neither the collimated appearance of 
the class outflows or the wide-angle morphology characteristic of the class I sources. Class 
II sources are generally found in regions of low ambient density (Myers et al. 1987), and the 
clumpy structures seen in Figure 1 may arise only in the (limited) regions where molecular 
gas around the nascent star has not yet been cleared by the protostellar outflow. RNO 91 
is a case in point. A large scale mosaic of 12 CO interferometer observations, covering an 
area about 5 x 10 4 by 10 5 AU, demonstrates that the outflow has produced a 3 x 10 4 AU 
wide cavity, with an opening angle of 160° (Lee et al. 2002; Lee & Ho 2005). Thus, our 
observations of the RNO 91 source, with a primary beam of ~ 5 x 10 3 AU, cannot fully map 
the cavity walls seen in the maps of Lee et al. (2002), instead we can expect to detect only 
remnant molecular material that has not been swept up out of the field of view of our maps. 
A similar scenario applies for T Tau (Momose et al. 1996) and perhaps for GK/GI Tau. 

In Figure 5 we plot the opening angle, 9, of each of our observed outflows as a function 
of the source age, together with published data of other spatially resolved molecular outflows 
from low-mass YSOs, within 500 pc of the Sun, where the outflow opening angle near the 
source could easily be measured (see Table 9). In cases where the value of 9 depends on 
velocity, we have adopted the largest opening angle. It is immediately clear that there is a 
tendency for the outflow opening angle to increase with the evolutionary stage of the outflow 
source. A fit to our data alone yields log(#/deg) = (1.1 ± 0.2) + (0.16 ± 0.4)log(t/yr), quite 
consistent with the fit to the combined data, log(6>/deg) = (0.7±0.2) + (0.26±0.4)log(t/yr). 

Although an increase of outflow opening angle with age has been inferred by other 
authors (e.g., Richer et al. 2000), it has never been clear if the trend was a real effect 
or due to the varying beam sizes in independent studies of individual sources observed 
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at lower resolution single telescopes. Indeed, high-resolution observations indicated that 
the opening angle of the B5-IRS1 outflow is widening at about 0.006° yr" 1 (Velusamy & 
Langer 1998), but there was no suggesting evidence to confirm that this was a general trend. 
It is true that interferometer mapping of a sample of molecular outflows from sources at 
different evolutionary stages concluded that outflow widths increase with age (Lee et al. 2002). 
However, we argue that the outflow opening angle is a better probe of the outflow's impact 
on the envelope, since it is measured close to the protostar, where the outflow-envelope 
interactions are strongest. Outflow width, by contrast, may be measured anywhere along 
the outflow length, including regions far away from the circumstellar envelope. Our results 
demonstrate unambiguously that outflow opening angles increase with source age, implying 
a progressively stronger outflow-envelope interaction. 



5.1.1. Why should we expect outflow widening? 

There is as yet no generally accepted explanation for the observed outflow widening. 
It has been proposed that jet axis wandering (precession) could produce wider cavities as 
the protostar evolves (Masson & Chernin 1993). However, although ample evidence for jet 
wandering exists (e.g., Reipurth et al. 1997; Terquem et al. 1999), in most sources the opening 
angle of the precession cone is smaller that the observed wide-angle outflow cavity (e.g., Yu 
et al. 1999; Reipurth et al. 2000; Paper I). Our study of the L 1228 outflow also shows that 
the detailed morphology of the wide-angle cavities is incompatible with their being produced 
by jet entrainment (Paper I). An attractive alternative is an underlying stellar wind with 
both a collimated and a wide-angle component (e.g., Kwan & Tademaru 1995; Konigl & 
Pudritz 2000, Shu et al. 2000; Matt et al. 2003) where the observed molecular outflow is 
predominantly driven by one of the two components depending on the age of the protostar. 
The X-wind model of Shang et al. (1998), for example, predicts a stellar wind with a high 
density along the axis (very similar to a jet) but decreasing density at increasing angles 
from the axis. It is conceivable that at very early ages only the dense collimated part of 
the wind can break out of the surrounding dense (infalling) envelope (e.g., Wilkin & Stahler 
2003), thereby producing the jet-like outflows observed in class sources. As the envelope 
loses mass, through infall onto the star-disk system and outflow entrainment along the axis 
(see § 5.4), the less dense and wider wind component will break through, entraining the gas 
unaffected by the collimated component. At this stage there will be almost no gas along the 
outflow axis, and the gas entrained by the wide-angle component of the wind will be the 
dominant structure in the molecular outflow, as seen in class I sources. A decrease in infall 
as the protostar evolves (e.g., Henriksen et al. 1997) could also allow the wind to carve out 
wide-angle cavities (see Delamarter et al. 2000). 
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5.2. Evolution in the Envelope's Morphology and Kinematics 

There also appears to be a change in the C 18 morphology as the outflow source evolves. 
From Tables 5 and 8 it is clear that for class sources the C 18 is usually elongated along 
an axis within 15° of that of the outflow. In Figure 6, two Class sources, RN043 and 
IRAS 3282, display velocity gradients along the (elongated) major axis of the envelopes, in 
the same sense as the outflow. There is no clear velocity gradient in the C 18 envelope of 
HH 114mms, although the gas at the extreme southwest edge of the envelope is blueshifted 
— similar to the CO outflow. Based on our detailed analysis of C 18 emission in RN043 
(Paper II), we suggest that the outflow is responsible for the velocity field and morphology 
of the high-density envelope around at least two class sources in our sample. 

On the other hand, Figure 6 and Tables 5 and 8 indicate the C 18 envelopes around the 
class I sources tend to be elongated perpendicular to the outflow axis, and there is no obvi- 
ously shared velocity distribution. For both RNO 129 and L 1228 the C 18 velocity gradient 
is aligned with the major axis of the envelope and perpendicular to the outflow axis, consis- 
tent with rotation. In HH300, the envelope velocity gradient is along the minor axis, parallel 
to the outflow axis, suggesting a flattened infalling envelope, where infall motions produce 
blueshifted and redshifted emission in the far and near sides of the structure, respectively 
(as seen in other sources, e.g., Torrefies et al. 1995; Ohashi et al. 1997; Momose et al. 1998). 
Contrary to what is observe in the class sources shown here and in the literature where the 
outflow is clearly impacting the dense circumstellar gas (i.e., Gueth et al. 1997; Wiseman et 
al. 2001; Beltran et al. 2004), the C 18 emission in HH300 is elongated perpendicular to the 
outflow at all velocities, and does not show any protrusion along the outflow axis. We are, 
therefore, certain that the HH300 velocity distribution is mainly (if not all) due to infall. 

For class II sources, either no C 18 emission is detected or the emission is offset in 
both position and velocity from the protostar. Specifically, the C 18 emission associated 
with RNO 91 is about 8" (1300 AU) southwest of the continuum peak and is seen only at 
velocities that are blueshifted with respect to that of the ambient envelope velocity. Similarly, 
the C 18 emission associated with T Tau is concentrated 15" (2100 AU) east of the young 
star and most of the emission emerges at velocities blueshifted by 0.2 to 0.5 km s _1 with 
respect to v env . In both cases the C 18 emission is concentrated toward the edge of the 
outflowing 12 CO structure seen in Figure 1. No C 18 emission within our sensitivity levels 
was detected from GK/GI Tau. 

As we describe in §5.4, we interpret these changes in the C 18 envelope properties 
between class and class II sources as a real evolutionary progression due to the protostellar 
outflow's effect on the dense circumstellar gas. 
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5.3. Envelope Mass Evolution 

In Figure 7 we plot our estimate of the envelope mass derived from the C 18 emission, 
■M-env 5 a function of source age, determined from the T^-age relationship of Ladd et 
al. (1998). Although the correlation is not very tight, M env is certainly greater for younger 
protostars than for more evolved YSOs. We would not expect any relation between T bo i and 
M env a priori since is obtained from the SED, and thus is an indirect measure of the 
dust content close to the protostar rather than the gas. 

A linear fit to the log(M em) )-log(age) plot results in: 

log(M em ,/M ) = (0.87 ± 0.5) - (0.45 ± 0.09)log(t/yr), (1) 

with a correlation coefficient of 0.8. Taking the derivative, leads to a rough estimate of the 
envelope mass- loss rate, M env , from: 

Menv ~ "3r L45 (2) 

Evidently, M env decreases rapidly with time from M env ~ 10~ 4 M yr -1 for sources with 
age 10 3 yrs, to 5 x 10~ 6 , and 2 x 10~ 7 M yr" 1 for ages of 10 4 , and 10 5 yrs, respectively. 
This is consistent with the core mass-loss rate found by Ladd et al. (1998), M oc t~ L52±0 - 55 , 
from lower resolution C 18 O(l-0) observations of a sample of Taurus cores harboring deeply 
embedded YSOs. 

Does the protostellar outflow play a role in the mass-loss of the circumstellar envelope? 
As discussed in §4.3 (and in Papers I and II in more detail), high- velocity 13 CO traces dense 
envelope gas that has been entrained by the protostellar wind and thus provide an estimate 
of the outflow rate of the relatively dense gas that could escape the envelope's gravitational 
potential. The dense gas mass outflow rate, M dense , shown in Table 6 is approximately con- 
stant across our sample, mostly ranging from only 1 to about 5 x 10~ 6 M yr -1 , independent 
of the source's age. Only two sources, IRAS 3282 and GK/GI Tau, deviate, with M dense 
values of ~ 4 and 2 x 10~ 7 M yr" 1 , respectively. The median (and average) value of M dense 
from our whole sample is between 1 and 2 x 10~ 6 M yr" 1 . This is very close to the value 
of the envelope mass-loss rate for late class 0/early class I sources with ages of about 10 4 
yrs, M env ~ 5 x 10~ 6 M yr -1 (see Equation 2). It is, therefore, conceivable that a fraction 
of the mass-loss in the envelope is due to the outflow entrainment of dense envelope gas. 
In class I (and older) sources, with ages > 10 5 yr, the envelope mass-loss rate is even less, 
M env < 2 x 10~ 7 M yr -1 , and is about a factor of ten smaller than the average M dense in our 
sample. Therefore, in class I (and older) sources even if only 10 to 20% of the dense mass 
entrained by the outflow eventually escapes from the envelope's gravitational potential, the 
outflow-envelope interaction can contribute significantly to the envelope's mass loss. 
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We note that in sources younger than ~ 10 3 yr the outflow may be responsible for only 
a small fraction of the total mass-loss in the envelope since the estimated dense gas outflow 
mass rate is less than 10% of the estimated envelope mass-loss rate. This accords well with 
theoretical and observational work suggesting that during this very early stage of protostellar 
evolution the envelope mass loss rate is dominated by gas infall onto the star-disk system 
(e.g., Ladd et al. 1998; Henriksen et al. 1997; Bontemps et al. 1996). 

5.4. Outflow-Envelope Interactions: An Empirical Model 

The widening of outflow cavities, the change in circumstellar envelope morphology and 
kinematics, and the decrease of circumstellar mass and mass-loss rate with source age all 
indicate that the protostellar outflow has a major role in the early stages of star formation, 
particularly in the evolution of the dense circumstellar envelope. Our data support the follow- 
ing empirical model of the evolution of outflow-envelope interactions (shown schematically 
in Figure 8). 

At very early stages of protostellar evolution, the young, powerful, and well collimated 
outflow entrains and pushes on gas as it burrows through the envelope, altering the distri- 
bution and kinematics of dense gas. This interaction causes the elongated distribution of 
the dense gas, as traced by the C 18 in our observations, and a velocity gradient along the 
outflow axis on scales of 10 3 AU (see Figure 8). At this early stage the envelope mass-loss 
rate is dominated by the infalling envelope gas onto the star-disk system. 

With time, the envelope mass decreases and the outflow cavities widen. This enables 
circumstellar gas at larger angles from the outflow axis to be entrained and cleared away 
by the outflow until, eventually, most of the dense circumstellar gas is concentrated outside 
the outflow cavities in structures perpendicular to the outflow axis, thereby decreasing the 
volume of the infalling region (see Figure 8). At this stage most of the remaining dense 
circumstellar gas has not yet been perturbed by outflow motions, and the velocity structure 
of the envelope is dominated by either infall or rotation, depending on the age of the source 
(see Hogerheijde 2001). The outflow-envelope interaction is mainly traced by the 13 CO 
outflow emission, as the lower density in the outer envelope precludes the use of C 18 as a 
probe. From this point (class I and later), the mass- loss rate of the envelope is dominated 
by the outflow's clearing effect. 

As the nascent star evolves further, the outflow lobes keep widening and by the time 
the source reaches class II status most of the dense circumstellar envelope has been swept 
away by the outflow or has fallen into the circumstellar disk-star system. Only remnants of 
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the original protostellar environment will be detectable (see Figure 8). In some cases, gas 
that has been entrained by the outflow may form relatively dense shells and clumps at the 
edge of the outflow lobe, that can be observed in 13 CO and C 18 emission around class II 
sources (as seen in our sample and in Welch et al. 2000). We speculate that as the YSO 
evolves further and reaches the end of the class II stage the very wide outflow cavities will 
finish clearing the surrounding gas that is left, halting the infall process altogether. 



5.5. Outflow Chemistry 

Our multi-molecular line observations also provide evidence that protostellar outflows 
can have a chemical impact on their envelope. As discussed in § 4.5 and shown in Table 7, 
we can estimate the HCO + abundance, Xhco+> m the outflow. Comparing our results 
with a "standard" abundance of ~ 7 x 10~ 9 (J0rgensen et al. 2004) we see that the HCO + 
abundance is clearly enhanced in the IRAS 3282 and HH114mms outflows. In the case of 
L1228, the lower limit of X HCO + in the blue and red outflow lobes is similar to, and less 
than a factor of two less than, the "standard" HCO + abundance, respectively. In the case of 
RNO 43, the estimated X HC o+ lower limit of the RNO 43 outflow is similar to the "average" 
HCO + abundance in envelopes around Class sources found by J0rgensen et al. (2004). It is 
possible that in these two sources the HCO + abundance may be slightly enhanced, but our 
data is not sufficient to firmly conclude this. In the other sources in our sample, the HCO + 
is similar to, or below, the standard abundance. 

The HCO + emission seen in Figure 3 arises mainly along the outflow cavity walls or 
the outflow lobes, and the observed overabundance is consistent with chemical models that 
predict HCO + enhancement in outflow walls due to the shock-triggered release of molecules 
from dust mantles and subsequent chemical reactions in the warm gas (Rawlings et al. 2000; 
Viti et al. 2002; Rawlings et al. 2004). Hence, it is clear that some of the sources in our 
sample power outflows that are affecting the chemical composition of their surroundings. On 
the other hand, HCO + outflows with abundances at or below the standard value must arise 
in dense (n > 10 3 ) envelope gas (see Evans 1999) that has been entrained by the outflow, 
where shocks have not seriously affected the chemical properties of the environment. In some 
sources the HCO + outflow emission does not trace the entire extent of the outflow wall (e.g., 
HH300, L1228) or does not trace both outflow lobes (e.g., HH114mms, RNO 43). This 
could result from differences in the environmental conditions such as density, shock-induced 
radiation, or chemical composition, or to the effects of optical depth. 
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6. Summary 

Our systematic, high angular resolution (3" to 7"), multi-molecular line, survey of the 
circumstellar environment within ~ 10 4 AU of nine protostars at different evolutionary 
stages has enabled detailed studies of the kinematics, morphology, density distribution, and 
chemistry of the circumstellar gas. These, in turn have permitted investigations of how the 
circumstellar envelopes and protostellar outflows change with time. 

Our 12 CO images trace high- velocity outflowing molecular material around all sources. 
We detect a clear trend in the morphology of the protostellar molecular outflows. The 
youngest sources (class 0) in our sample power molecular outflows with jet-like morphologies 
or cone-shaped lobes with opening angles less than 55°. The sightly more evolved (class I) 
sources in our sample have molecular outflows with lobe opening angles of more than 75°. 
Outflows from the most evolved young stars in our sample, the class II sources, have even 
wider lobes or no definite shape or structure. Combining our data with that from a number 
of sources in the literature shows that outflow opening angle close to the source widens with 
time. 

The denser, lower velocity gas probed by our 13 CO and HCO + observations appears to 
arise from the outflow cavity walls or very close to the protostar. We suggest that this 
is dense gas from the outer regions of the circumstellar envelope that has been entrained 
by the high velocity flow, thereby eroding the envelope and helping to widen the outflow 
cavities. We also suggest that evolutionary changes in the morphology and velocity field 
of the dense circumstelar envelopes, detected in our C 18 images, are mainly caused by 
outflow-envelope interactions. Class source envelopes, for example, are elongated with 
velocity gradients along the outflow axis, indicating that the young and powerful protostellar 
outflows can entrain dense envelope gas. Class I envelopes by contrast are elongated more or 
less perpendicular to the outflow axis and concentrated outside the outflow lobes, as if most 
of the dense gas along the outflow axis has already been displaced. By the class II stage 
the sparse C 18 emission observed is constrained to the outflow lobe edges, consistent with 
most, if not all, of the dense gas being cleared. 

Not unexpectedly, our results show a decrease in envelope mass as the age of the pro- 
tostar increases. Comparisons of the estimated envelope mass-loss rate with the dense gas 
outflow rate imply that the protostellar outflow plays an important role in envelope mass- 
loss during and after the class I stage. For younger protostars, envelope mass-loss is more 
likely to be dominated by large infall rates. Finally, enhanced abundance of HCO + in the 
outflow in some of the sources in our sample indicate that the chemical composition of the 
environment around these protostars is affected by shock-induced chemical processes in the 
outflow-envelope interface. 
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Taken together, these trends demonstrate that outflows from protostars have a major 
physical and chemical impact on their environments throughout the entire star formation 
process and are an integral part of that process. An understanding of their origins, evolution, 
and detailed characteristics is clearly critical for theories of star formation. 

H.G.A. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship 
under award AST-0401568. The Owens Valley Radio Observatory millimeter array is sup- 
ported by the NSF grant AST-0228955. We are greateful to Chin-Fei Lee, Mario Tafalla and 
Jennifer Hatchell for sharing their single-dish data with us. 

A. Description of Sources 
A.l. HH114mms 

The HH 114mms protostar is a Class source in L1589, a cloud in the western part 
of the A Ori molecular shell. It was first detected at 1.3 mm (Chini et al. 1997), and 
has also been observed at 3.6 cm (Rodriguez & Reipurth 1996). The 12 CO(1-0) map of 
the nearby HH 114/115 molecular outflow by Lee et al. (2002) partially overlaps the blue 
lobe of the HH114mms CO outflow. Our data show the HH114mms CO outflow is highly 
collimated, with a small bend in both outflow lobes about 7" away from the source, where 
the outflow axis appears to change direction. This morphology may indicate that the axis of 
the underlying protostellar wind has changed over time. CO channel maps show that in both 
outflow lobes the highest velocity gas is concentrated farthest from the source. In addition, 
the opening angle of the CO blue lobe widens slightly at lower outflow velocities, rather like 
the HH211 outflow (Gueth & Guilloteau 1999). The morphology and velocity distribution 
of the outflowing gas suggests that an underlying jet-like wind drives the molecular outflow. 



A. 2. RN0 43 

A detailed study of the RNO 43 molecular outflow, the circumstellar envelope and the 
outflow-envelope interactions within 2 x 10 4 AU of the protostar is presented in Paper II. 
This outflow lies in the L 1582B dark cloud at a distance of about 460 pc, in the A Ori 
ring (Bence et al. 1996). The outflow arises from IRAS 05295+1247, a deeply embedded 
class source that has been observed from far-infrared to centimeter wavelengths (Cohen 
et al. 1985; Beichman et al. 1986; Zinnecker et al. 1992; Reipurth et al. 1993; Anglada et 
al. 1992; Chini et al. 1997). The 3.4 pc long RNO 43 optical outflow includes four major 
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Herbig-Haro (HH) knots, HH 243, HH 244, HH 245, and HH 179 (Cohen 1980; Reipurth 
et al. 1997). The large-scale molecular (CO) outflow extends about half a degree (~ 4 pc) 
on the sky (Cabrit et al. 1988; Bence et al. 1996). The positions of the HH knots and the 
morphology of the giant molecular outflow indicate that the axis of the RNO 43 flow has 
changed over time (Bence et al. 1996; Reipurth et al. 1997). 

A.3. IRAS 3282 

Single dish 12 CO(2-l) observations of the environment surrounding IRAS 3282 by Bachiller 
et al. (1991) revealed a young, 0.5 pc long, fast, and highly-collimated bipolar molecular out- 
flow. The high velocity outflow (~ 50 km s _1 ) gas has a jet-like morphology, while the slower 
(~ 10 km s _1 ) component exhibits wider outflow lobes. Observations of the ammonia core 
suggest that the dense core gas is being accelerated along the outflow cavity walls (Tafalla 
et al. 1993). High resolution OVRO observations of the 1.3 mm dust continuum, reveal that 
IRAS 3282 is a binary source with a 420 AU separation and a mass ratio of 0.2 (Launhardt 
2004). Our 12 CO(1-0) data show that at outflow velocities between about 2 and 7 km s _1 
both lobes have a V-shaped morphology, consistent with a jet bow shock entrainment model 
(e.g., Gueth et al. 1996; Gueth & Guilloteau 1999). 

A.4. HH300 

IRAS 04239+2436, the source of the HH300 outflow, is located in the westernmost re- 
gion of the B18 cloud (sometimes refereed as the B18w region) in the Taurus molecular cloud 
complex. Infrared photometry and spectra indicate that the source is a deeply embedded 
class I protostar (Myers et al. 1987; Greene & Lada 1996). Dust continuum emission sur- 
rounding the source was mapped at 1.3 mm (at ~ 12" resolution) by Motte & Andre (2001). 
The HH300 optical outflow is about 1.2 pc long and comprises three different HH knots 
(Reipurth et al. 1997). The density, velocity and momentum distribution of the large-scale 
molecular outflow indicates that it is driven by an underlying precessing and episodic wind 
(Arce & Goodman 2001). HST images at 1.644 and 2.122 /im by Reipurth et al. (2000) 
indicate that IRAS 04239+2436 is a binary source. In addition, they reveal the existence 
of a parabolic nebula with a small (1000 AU) jet along its symmetry axis, northeast of the 
source (i.e., the blueshifted side of the outflow). In our 12 CO images, the blue and red molec- 
ular outflow lobes have shapes and orientations similar to that of the near-infrared parabolic 
nebula. 
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A.5. L1228 

A detailed study of the infalling envelope, the wide-angle L 1228 molecular outflow, 
and its impact on the circumstellar envelope is presented in Paper I. This outflow is named 
after its parent molecular cloud, located in the Cepheus flare cloud complex at a distance 
of 200 pc (Kun 1998). The outflow source, IRAS 20582+7724, is an an embedded (Class 
I) low-luminosity protostar associated with a ~ 1 pc-long bipolar molecular outflow, first 
observed in CO(l-O) by Haikala & Laureijs (1989). IRAS 20582+7724 also powers a ~ 1.5 pc- 
long optical flow, HH 199 (Bally et al. 1995). The single-dish, multi-molecular line study by 
Tafalla & Myers (1997) indicates that the dense core gas is being disrupted by the L1228 
outflow. The tentative detection of a binary companion is reported in Paper I. 

A.6. RN0129 

RNO 129 is located about 1° north of the L 1228 core and is believed to be associated 
with the same cloud (Kun 1998). The protostellar source associated with the RNO 129 
nebulosity is IRAS 21004+7811 (Ogura & Sato 1990). An extensive recent study of the 
region by Movsessian & Magakian (2004) presents optical images of the RNO 129 nebulosity, 
and the HH objects near it, and concludes that the central star is a binary that powers a 
collimated (optical) outflow. Here we show the first maps of the molecular outflow associated 
with RNO 129. 



A. 7. RNO 91 

RNO 91 and its source, IRAS 16316-1540, are located in the L43 cloud, part of the 
p Oph molecular cloud complex at a distance of about 160 pc. Optical, CO, CS, and NH3 
observations of the core suggest that the molecular outflow has blown through the dense 
gas, creating a cavity south of the outflow source (Mathieu et al. 1988). Based on single- 
dish and interferometer maps at higher angular resolution, the molecular outflow has been 
successfully modeled as a slowly expanding shell driven by a wide-angle wind (Bence et 
al. 1998; Lee et al. 2000, 2001; Lee & Ho 2005). Given the large width of the outflow near 
its base, about 10 4 AU, our observations, which are limited to regions within about 6.5 x 10 3 
AU of the source, do not show the full impact of the RNO 91 outflow on the surrounding 
envelope. Fortunately, the mosaic of interferometer observations by Lee & Ho (2005), which 
cover the surrounding environment within at least 10 4 AU of IRAS 16316-1540, show these 
outflow-envelope interactions. 
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A.8. T Tau 

The T Tau binary system comprises the optically visible star T Tau N and an infrared 
companion, T Tau S, separated by only 0.7". Both binary elements are known outflow 
sources. The outflow mostly likely driven by T Tau S has an approximately north-west axis 
(e.g., Solf & Bohm 1999; Reipurth et al. 1997, and references therein), while the other seen 
nearly pole-on, with a northeast-southwest or east-west axis, is most likely driven by T Tau 
N (Solf & Bohm 1999; Momose et al. 1996, and references therein). It is generally accepted 
that both T Tau N and T Tau S are of similar age; T Tau N is optically visible because its 
light can escape through the nearly pole-on outflow cavity walls (e.g., Momose et al. 1996), 
while T Tau S is obscured by an envelope or the outer parts of T Tau N's disk, and is highly 
reddened (e.g., Hogerheijde et al. 1997b). However, the bolometric temperatures of the two 
stars are extremely different: 501 and 3452 K for T Tau S and T Tau N, respectively (Chen 
et al. 1995). Since we do not have separate bolometric temperature measurements for each 
element of other close (< 5") binary systems in our sample (i.e., IRAS 3282, HH300, L 1228, 
RN0129), and since T Tau S and T Tau N should be coeval, we adopt a single T bol for 
both stars using flux measurements at different wavelengths that do not resolve the binary 
companions. 

The molecular gas morphology and kinematics in the region surrounding T Tau are 
quite complex, and different observational studies offer different interpretations. Weintraub 
et al. (1989b) conducted the first interferometer observations (with a beam of ~ 6") of the 
12 CO(1-0) and 13 CO(1-0) surrounding T Tau and interpreted their observations to show a 
circumstellar Keplerian disk. The combined interferometer and single-dish 13 CO(1-0) ob- 
servations of Momose et al. (1996) show a different picture, and were interpreted to show 
biconical outflowing shells. The interpretation by Schuster et al. (1997) based on their single- 
dish 12 CO(3-2) and C 18 0(2-l) observations (with beams of 12 to 14") is similar to that of 
Momose et al. (1996), but they claim more than one outflow is needed to explain the ob- 
served gas morphology and kinematics. Observations of the extended molecular hydrogen 
in the v = 1 — S(l) line show the emission appears to arise from a precessing jet in an 
approximate north-south axis powered by T Tau S (van Langevelde et al. 1994b). Interfer- 
ometer observations of the HCO + (1-0) towards T Tau show possible evidence of a infalling 
circumstellar envelope (van Langevelde et al. 1994a). 

The OVRO 12 CO(1-0) maps presented here are much more sensitive, have a better 
angular resolution (by a factor of two), and better uv-coverage (by a factor more than 4) 
than those presented by Weintraub et al. (1989b). Our data does not show the Keplerian 
circumstellar disk suggested by Weintraub et al. (1989b). Instead, at velocities far away from 
the cloud velocity (i.e., V[ sr < 6.5 km s _1 and Vi sr > 9.5 km s _1 ) the 12 CO emission shows 
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outflow lobes with an east-west or northeast-southwest axis, consistent with the Momose et 
al. (1996) picture — where most of the outflowing molecular emission comes from an outflow 
powered by T Tau N. Close to the cloud velocity the 12 CO emission (not shown here) is 
complex and extended, and it is clear that we are missing a substantial amount of flux. 
Our 13 CO, C 18 and HCO + maps are consistent with the lower angular resolution maps 
presented by Momose et al. (1996), Schuster et al. (1997), van Langevelde et al. (1994a), 
respectively. 

A.9. GK/GI Tau 

GK Tau and GI Tau are two classical T Tauri stars located in the B18 cloud in Taurus, 
that are only separated by 13" (~ 1820 AU). In addition, GK Tau has a faint companion 
only 2.5" (350 AU) away (Reipurth & Zinnecker 1993). Estimates of the age and bolometric 
temperature for GI Tau and GK Tau indicate that these stars are coeval (Hartigan et al. 1994; 
Chen et al. 1995). Here we assume a common bolometric temperature of 2700 K. Aspin & 
Reipurth (2000) detected several HH knots near the triple system, and were able to determine 
that GK Tau (or its faint companion) is the driving source of at least some of the HH knots. 
Their observations also reveal that the entire triple system is surrounded by a common 
envelope of reflection nebulosity. Our CO maps of the region are the first to show the 
molecular outflow associated with the GK/GI system. 
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Fig. 1. — Gallery of 12 CO(1-0) outflows. Blue (red) contours show the blueshifted (red- 
shifted) emission from the outflow lobes of the sources in our sample, integrated over the 
velocity ranges given in Table 4. The upper, middle and lower panels show the class 0, class 
I, and class II objects, respectively. A cross marks the position of the protostar given by the 
millimeter continuum emission peak. In each map, the synthesized beam is in the lower left 
corner. A scale equivalent to 10000 AU at the assumed distance of each source is also given. 
The arrowed dashed line in each panel represents the presumed outflow axis. The values of 
the first contour and the subsequent contour steps for the blue (red) lobe of each outflow are 
given inside brackets in units of Jy km s^ 1 , at the lower left (right) corner of each panel. 

Fig. 2. — Gallery of 13 CO(1-0) outflows. Solid blue (red) contours show the blueshifted 
(redshifted) integrated intensity. The dashed contours show the 12 CO outflows (shown in 
Figure 1). The velocity range of integration for each 13 CO outflow lobe is given in Table 4. 
The cross symbol indicates the position of the protostar, given by the millimeter continuum 
emission peak. The synthesized beam of the 13 CO ( 12 CO) map is shown on the lower right 
(left) corner of each panel. The scale in each panel shows the equivalent of 5000 AU at the 
assumed source distance. The values of the first contour and the subsequent contour steps 
for each outflow lobe are given inside brackets in units of Jy km s" 1 . In the case of RN0 43, 
the "blueshifted" lobe is shown in grey, as the gas velocities are very close the ambient 
envelope velocity. The RN043 outflow axis is very close to the plane of the sky, thus this 
13 CO structure most likely traces the walls of the outflow cavity of the blueshifted outflow 
lobe, albeit its velocities (see Paper II). 

Fig. 3. — Gallery of HCO + (1-0) outflows. Solid blue (red) contours represent integrated 
intensity of blueshifted (redshifted) HCO + emission. Dashed contours represent 12 CO molec- 
ular outflow. The synthesized beam of the HCO + ( 12 CO) map is shown on the lower right 
(left) corner of each panel. The values of the first contour and the subsequent contour steps 
for the blue (red) lobe of each outflow are given inside brackets in units of Jy km s _1 , at the 
lower left (right) corner of each panel. 

Fig. 4. — C 18 integrated intensity (solid green) contour maps superimposed on 12 CO out- 
flows (dashed contours). The velocity range of integration for each source is given in Table 8. 
The synthesized beam of the C 18 ( 12 CO) map is shown on the lower right (left) corner of 
each panel. The scale in each panel shows the equivalent of 5000 AU at the assumed source 
distance. The values of the first contour and the subsequent contour steps are given in 
brackets at the bottom of each panel. 

Fig. 5. — Outflow opening angle as a function of source age. Source age was obtained using 
the T bo ;-age relation from Ladd et al. (1998). Filled circles show sources in our sample for 
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which we were able to measure an outflow opening angle (see Tables 1 and 5), and the filled 
triangles are sources from the literature (see Table 9). The dashed line represents the fit to 
all the data shown in the figure. The errors in the opening angle come from the uncertainty 
of the angle measurement, while the errors in the age come mainly from the T^-age relation. 

Fig. 6. — First moment maps of the C 18 O(l-0) line from our class (top row) and class 
I (bottom row) sources. The first moment between two velocities, v\ and i>2, is defined as 
J^TbvcIv/ f^Tgdv, where Tg is the brightness temperature. The color bar to the right 
of each panel shows the color that represents the velocity in each map. Contours represent 
integrated intensity levels similar to those of Figure 4. The arrow shows the 12 CO outflow 
axis, from Figure 1. 

Fig. 7. — Envelope mass as a function of protostar age. The envelope mass was obtained 
from the C 18 emission (see Table 8). The protostar age was obtained using the T^-age 
relation from Ladd et al. (1998). The dashed line represents the fit to all the data shown 
in the figure. The error bars in the envelope mass reflect the 30% calibration uncertainty in 
the data, while the error in the age come mainly from the Tfe ;-age relation. 

Fig. 8. — Schematic picture of outflow-envelope interaction evolution. Dark grey regions 
denote high-density (envelope) gas, mostly traced by C 18 in our observations. Light grey 
regions show the molecular outflow traced by the 12 CO and 13 CO. Arrows indicate the gas 
motion. See 5 5.4 for details. 
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Table 1. Observed Sources 



Source Name 


Alternate Name 


Observation 


Coordinates 


d 


Class 


T boi a 


SED References' 3 






R.A. [2000] 


Dec. [2000] 


[PC] 


[K] 






HH 1 14mms 




5 A 18 m 15 s .4 


7°12'00" 


460 





40 


1 


IRAS 05295+1247 


RNO 43 


5 32 19.4 


12 49 42 


460 





40 


2,3 


IRAS 03282+3035 


IRAS 3282 


3 31 21.0 


30 45 31 


350 





50 


2,4,5 


IRAS 04239+2436 


HH 300 


4 26 56.3 


24 43 36 


140 


I 


225 


2,5,6,7 


IRAS 20582+7724 


L1228 


20 57 13.0 


77 35 47 


200 


I 


345 


2,7,8 


IRAS 21004+7811 


RNO 129 


20 59 13.8 


78 23 05 


200 


I 


650 


2,7 


IRAS 16316-1540 


RNO 91 


16 34 29.3 


-15 47 01 


160 


II 


715 c 


2,7,9 


T Tau 




4 21 59.4 


19 32 07 


140 


II 


1235 


2,7,9,10,11,12,13 


GK/GI Tau d 




4 33 34.4 


24 21 07 


140 


II 


2700 c 


7,9 



a Bolometric temperature calculated from spectral energy distribution following Chen et al. (1995) 
and combining our new observations with published data. Note that for T Tau, the flux measure- 
ments used do not resolve the individual components (see A. 8). 

b References used to obtain spectral energy distribution for each source: (1) Chini et al. (1997); 
(2) IRAS Catalog of Point Sources; (3) Zinnecker et al. (1992); (4) Dent et al. 1998; (5) Motte 
& Andre (2001); (6) Myers et al. (1987); (7) 2MASS All-Sky Catalog of Point Sources; (8) Arce 
& Sargent (2004); (9) Herbig & Bell (1988); (10) Weintraub et al. (1987); (11) Weintraub et 
al. (1989a); (12) Beckwith et al. (1990); (13) Beckwith k Sargent (1991) . 

c Bolometric temperature used here was taken directly from Chen et al. (1995). 

d GK Tau and GI Tau are less than 15" apart and both are inside our field of view. 
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Table 2. Characteristics of Line and Continuum Maps 



Source Line Beam Size Map RMS a Line/Continuum Beam Size Map RMS a 

[arcsec] [Jy beam -1 ] [arcsec] [Jy beam -1 ] 

HH 114mms 



RNO 43 



IRAS 3282 



HH300 



L1228 



RNO 129 



RNO 91 



L Lau 



12 co 


1-0) 


4.3 x 3.9 


0.08 


HCO+(1-0) 


6.4 x 5.4 


0.11 


13 co 


1-0) 


5.1 x 5.0 


0.09 


2.7 mm cont. 


6.4 x 5.9 


0.003 


L (J 


l-UJ 


5.3 X 4.9 


U.U7 


3.4 mm cont. 


b.3 x 6.1 


U.UUU7 


12 co 


1-0) 


4.4 x 4.0 


0.10 


HCO+(1-0) 


5.4 x 5.2 


0.13 


13 co 


1-0) 


5.1 x 4.1 


0.10 


2.7 mm cont. 


5.1 x 5.0 


0.0007 


rU8n 
L (J 


l-UJ 


A *7 w A 1 

4.7 X 4.1 


n no 

U.U8 


3.4 mm cont. 


5.6 X 4.6 


n nnn A 

U.UUU4 


12 co 


1-0) 


3.5 x 2.9 


0.06 


HCO+(1-0) 


6.4 x 5.9 


0.13 


13 co 


1-0) 


4.9 x 4.1 


0.08 


2.7 mm cont. 


4.3 x 4.0 


0.001 


u u 


1-U J 


o.o X o.4 


U.Ub 


3.4 mm cont. 


O.U X 4.U 


U.UUUb 


12 CO 


1-0) 


4.4 x 3.2 


0.07 


HCO+(1-0) 


6.9 x 6.6 


0.11 


13 CO 


1-0) 


4.6 x 3.3 


0.10 


2.7 mm cont. 


4.8 x 3.7 


0.0008 


rU8n 


'i n\ 
1-U J 


4.1 X 3.0 


U.Ub 


3.4 mm cont. 


4.8 X 3.8 


n nnn A 
U.UUU4 


^CO 


1-0) 


4.3 x 3.8 


0.08 


HCO + (1-0) 


7.1 x 6.6 


0.11 


13 co 


1-0) 


4.8 x 4.0 


0.07 


2.7 mm cont. 


6.0 x 5.7 


0.0006 


c 18 o 


1-0) 


4.7 x 3.9 


0.08 


3.4 mm cont. 


6.3 x 5.5 


0.0004 


12 co 


1-0) 


5.7 x 4.4 


0.11 


HCO+(1-0) 


8.5 x 5.3 


0.10 


13 co 


1-0) 


6.4 x 5.0 


0.10 


2.7 mm cont. 


7.9 x 6.9 


0.002 


c 18 o 


1-0) 


6.8 x 5.7 


0.08 


3.4 mm cont. 


9.1 x 6.0 


0.0007 


12 co 


1-0) 


6.2 x 4.5 


0.16 


HCO+(1-0) 


11.7 x 6.3 


0.25 


13 co 


1-0) 


6.1 x 4.9 


0.15 


2.7 mm cont. 


6.6 x 4.7 


0.001 


c 18 o 


1-0) 


6.6 x 5.0 


0.12 


3.4 mm cont. 


13.3 x 7.8 


0.0009 


12 co 


1-0) 


3.4 x 3.0 


0.08 


HCO+(1-0) 


5.2 x 4.5 


0.12 


13 co 


1-0) 


3.7 x 3.6 


0.09 


2.7 mm cont. 


4.6 x 4.3 


0.002 


c 18 o 


1-0) 


4.3 x 3.8 


0.06 


3.4 mm cont. 


4.9 x 4.1 


0.0005 
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Table 2 — Continued 



Source Line Beam Size Map RMS a Line/Continuum Beam Size Map RMS a 

[arcsec] [Jy beam -1 ] [arcsec] [Jy beam -1 ] 



GK Tau 














12 CO(1-0) 


5.9 x 4.4 


0.10 HCO+(1-0) 


5.3 x 4.1 


0.16 






13 CO(1-0) 


7.0 x 6.1 


0.12 2.7 mm cont. 


7 1 v K 3 


n nnns 






C 18 O(l-0) 


6.9 x 5.7 


0.09 3.4 mm cont. 


0.4 x 6.1 


0.0U04 






a In line maps, the rms per velocity channel is shown. 












Table 3 


. Continuum Results 










Source R.A. 


Dec. 


Z? a 
-^2.7 


T? a 
-f 3.4 




T d c 


M cont d 






[mJy] 


[mJy] 




[K] 


[M ] 


HHlWmms 5' l 18 m 15 s .2 


7°12'02" 


110 ±20 


57 ±11 


0.4 ±0.2 


30 


0.24 


RN0 43 5 32 19.4 


12 49 41 


26 ±6 


11 ±2 


1.0 ±0.1 


35 


0.19 


IRAS 3282 3 31 20.9 


30 45 30 


61 ± 12 


30 ±6 


0.76 ±0.2 


28 


0.17 


HH300 4 26 56.3 


24 43 35 


14 ±2 


8±2 


1.5 ±0.4 


40 


0.02 


L1228 20 57 12.9 


77 35 44 


45 ±9 


21 ±4 


1.2 ±0.4 


33 


0.10 


RNO 129 20 59 13.6 


78 23 03 


62 ± 13 


28 ±6 


1.0 


40 


0.07 


RN0 91 16 34 29.3 


-15 47 01 


8±2 


3 ±0.5 


1.0 


50 


0.005 


T Tau 4 21 59.5 


19 32 06 


60 ± 12 


33 ±6 


0.25 ±0.2 


45 


0.005 


GK/GI Tau 4 33 34.53/34.13 24 21 06/18 < 2.0 f /4.0 ± 0.8 1.0 ± 0.2/3.0 ± 0.6 


1.0 


40 


0.001/0.00 



a Measured continuum fluxes at 2.7 mm (-F2.7), and at 3.4 mm (-F3.4). 

b We show the 1-a error of our estimate of (5, the power-law index of the dust emissivity. We assigned 
(3 = 1 to sources for which we do not have enough data to estimate its value (see text for details). 

Uncertainties in Td are about 20%. 

Uncertainties in M cont are about 30%. 

e The dust continuum emission associated with GI Tau (a companion star to GK Tau) is also detected in 
the GK Tau field 

'Quoted value is 3-cr upper limit. 
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Table 4. Velocity Ranges of Molecular Outflow Lobes 



Source 12 CO Outflow 13 CO Outflow HCO+ Outflow 

Blue Lobe Red Lobe Blue Lobe Red Lobe Blue Lobe Red Lobe 
km s _1 km s _1 km s -1 km s _1 km s _1 km s _1 



HHlWmms 


[-10.05, -3.87] 


[0.36, 8.16] 


[-1.26, -1.08] 


[0.10, 0.62] 




[0.19, 1.23] 


RN0 43 


[6.02, 8.95] 


[11.22, 16.10] 


[10.11, 10.45] a 


[10.63, 11.14] 




[10.76, 11.17] 


IRAS 3282 


[-2.75, 5.05] 


[8.63, 16.10] 


[5.98, 6.49] 


[7.18, 7.52] 


[4.28, 6.05] 


[7.11, 9.52] 


HH300 


[3.45, 6.05] 


[7.03, 10.28] 


[5.37, 6.02] 


[7.04, 7.90] 


[5.86, 6.07] 




L1228 


[-17.50, -9.05] 


[-6.13, 2.00] 


[-9.63, -8.60] 


[-6.90, -6.40] 


[-10.90, -9.02] 


[-6.92, -5.02] 


RNO 129 


[-17.75, -10.93] 


[-6.38, -2.14] 


[-9.87, -8.17] 


[-6.63, -5.79] 


[-8.84, -8.21] 


[-6.63,-5.90] 


RN0 91 


[-6.65, -0.80] 


[1.15, 2.13] 


[-0.70, -0.02] 


[1.01, 1.34] 


[-0.34, -0.10] 




T Tau 


[1.17, 6.70] 


[9.30, 13.53] 


[5.44, 7.14] 


[8.52, 10.05] 


[5.02, 7.19] 


[8.81, 9.81] 


GK/GI Tau 


[4.95, 6.25] 


[7.55, 9.50] 


[5.87, 6.38] 


[7.42, 7.92] 







a Although the entire range of velocities is not blucshifted with respect to the RNO 43 ambient envelope 
velocity, 10.2 km s , we are convinced the 13 CO structure at these velocities traces the outflow cavity walls 
of the RNO 43 blueshifted outflow (see Paper II). 



Table 5. Molecular Outflow Parameters I 



Source 


Outflow Axis PA 


Opening Angle a 


M out 


P b 

r out 


E ut h 




(deg. E of N) 


[deg.] 


(10- 3 M ) 


(10" 3 M Q km s- 1 ) 


(10 40 erg) 


HH 114mms 


60 


25 


30 


91 


484 


RNO 43 


60 


45 


150 


119 


142 


IRAS 3282 


125 


50 


15 


44 


195 


HH300 


59 


80 


8 


6 


2 


L1228 


79 


95 


13 


26 


92 


RNO 129 


95 


125 


24 


36 


94 


RNO 91 


155 


160 


11 


11 


20 


T Tau 


75 


100 


21 


25 


37 


GK/GI Tau 


c 




6 


4 


4 



a Outflow opening angle was estimated the following way: in HH114mms we used the lowest velocity 
blueshifted 12 CO; in RNO 43 and L 1228 we used the blushifted 13 CO outflow lobe shown in Figure 2; in 
IRAS 3282, HH300, RNO 129, and TTau we used both lobes of the 12 CO outflow shown in Figure 1; in 
RNO 91 we used the opening angle estimate from Lee & Ho (2005); in GK/GI Tau we unable to determine 
the opening angle from our maps. 

b Estimated value not corrected for outflow inclination with respect to the plane of the sky 

c Unknown outflow axis due to uncertainty as to which source powers the observed outflow. 
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Table 6. Molecular Outflow Parameters II 



Source 


M out a 


1 CO 


A/f a 
1V1 dense 




(1(T 7 Mq yr- 1 ) 


(10- 6 Mq km s- 1 yr" 1 ) 


(1(T 7 M yr- 1 ) 


HHlWmms 


41 


16 


11 


RN0 43 


73 


9 


51 


IRAS 3282 


31 


11 


4 


HH300 


18 


1 


16 


L1228 


29 


15 


13 


RNO 129 


56 


14 


37 


RN0 91 


16 


3 


11 


T Tau 


47 


8 


36 


GK/GI Tau 


2 


0.2 


2 



a Estimated value not corrected for outflow inclination with respect to the plane of the sky 



Table 7. HCO+ Abundance in Outflows 



Source [HCO+]/[H 2 ] 

Blue Lobe Red Lobe 



(lO" 9 ) (10- 9 ) 



HH 114mms 




> 20 


RNO 43 




> 2 


IRAS 3282 


> 2 


> 60 


HH300 


5 


1 


L1228 


> 7 


> 4 


RNO 129 


0.5 


0.6 


RNO 91 


1 




T Tau 


4 


3 


GK/GI Tau 
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Table 8. C 18 O(l-0) Envelope Parameters 



Source 


Velocity Range a 


Size b 


P.A. C 




Venv 


Reference d 




(km s _1 ) 


(AU) 


(deg. E of N) 


(M e ) 


(km s _1 ) 




HH 1 14mms 


[-1.09,0.28] 


2700 x 1400 


69 


0.13 


-0.46 


1 


RNO 43 


[9.26,11.48] 


6200 x 5000 


75 


0.87 


10.20 


1 


IRAS 3282 


[6.32,7.68] 


1900 x 1600 


113 


0.08 


6.96 


1 


HH300 


[5.84,7.21] 


3700 x 1300 


161 


0.03 


6.50 


2 


L1228 


[-8.77,-7.41] 


2400 x 1500 


161 


0.05 


-8.00 


1 


RNO 129 


[-8.35,-6.29] 


2000 x 1500 


149 


0.04 


-7.60 


1 


RNO 91 


[-0.35,0.33] 






0.01 


0.57 


1 


T Tau 


[7.15,8.34] 






0.02 


7.90 


3 


GK/GI Tau 








< 0.001 c 


6.90 


4 



a Lower and upper limits of velocity range for which C 18 emission was detected above the 3-cr 
level. 

b FWHM size of C 18 envelope from a gaussian fit to the integrated intensity map, deconvolved 
by the beam. 

c Deconvolved position angle of the gaussian fitted to the C 18 integrated intensity emission. 

d References for estimate of v env : (1) our observations; (2) Fuller & Ladd (2002); (3) average of 
values used by Momose et al. (1996) and van Langevelde et al. (1994a); (4) peak velocity at position 
of GK Tau from Mizuno et al. (1995). 

c Quoted value is 3-cr upper limit. 
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Table 9. Outflow Opening Angle from Literature 



Source 


(K) 


Reference 3, 


Opening Angle b Reference 
(deg) 


IRAM 04191 


18 


1 


45 


1 


HH211 


30 


1 


15 


2 


L1448 


56 


1 


40 


3 


L1527 


59 


1 


70 


4 


L1157 


62 


1 


35 


5 


HH 114/115 


65 


2,3,4 


30 


1 


RNO 40 


65 


2,3,4 


30 


6 


HH111 


70 


2,3,4 


20 


6 


L 1551-IRS5 


97 


1 


100 


7 


B 5-IRS1 


185 


4,5,6 


90 


8 



a Tfe ; obtained from estimates shown in (1) Buckle & Fuller (2002), or calculated from SED using 
data from the literature: (2) Reipurth et al. (1993); (3) Dent et al. (1998); (4) IRAS Catalog of 
Point Sources; (5) 2MASS All-Sky Catalog of Point Sources; (6) Langer et al. (1996). 

b Uncertainty in opening angle estimate is about 5 to 10°. 

c Reference for opening angle estimate: (1) Lee et al. (2002); (2) Gueth & Guilloteau (1999); (3) 
Bachiller et al. (1995); (4) Ohashi et al. (1997); (5) Gueth et al. (1996); (6) Lee et al. (2000); (7) 
Ohashi et al. (1996); (8) Langer et al. (1996). 
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